hep-ph/9811354 



FERMILAB-Pub-98 /366-T 
NHCU-HEP-97-02 



Vector Quark Model and B — ► X s "f Decay 



Chia-Hung V. Chang^ 1,2 ) Darwin Chang^ 1 ' 3 ), and Wai-Yee Keung^ 4 ' 3 ) 

W NCTS and Physics Department, National Tsing-Hua University, 
Hsinchu, 30043, Taiwan, R.O.C. 

^Physics Department, National Taiwan Normal University, Taipei, Taiwan, R.O.C. 

^Fermilab, P.O. Box 500, Batavia. Illinois, 60510, USA. 

^Physics Department, University of Illinois at Chicago, IL 60607-7059, USA 

Abstract 

We study the B meson radiative decay B — > X s ^ in the vector quark model. 
Deviation from the Standard Model arises from the non-unitarity of the charged current 
KM matrix and related new FCNC interactions. We establish the relation between the 
non-unitarity of charged current mixing matrix and the mixing among the vector quark 
and the ordinary quarks. We also make explicitly the close connection between this non- 
unitarity and the flavor changing neutral currents. The complete calculation including 
leading logarithmic QCD correction is carefully carried out. Using the most updated 
data and the NLO theoretical calculation, the branching fraction of the observed B 
meson radiative decay places a limit on the mixing angles as stringent as that from the 
process B — > X/j,£l. 



PACS numbers: 13.38.Dg, 12.15.Ff, 12.90. +b 



1 Introduction 



A simple extension of Standard Model (SM) is to enlarge the particle content by adding 
vector quarks, whose right-handed and left-handed components transform in the same way 
under the weak £77(2) x U(l) gauge group. This extension is acceptable because the anoma- 
lies generated by the vector quarks cancel automatically and vector quarks can be heavy 
naturally. Vector quarks also arise in some Grand Unification Theory (GUT). For example, 
in some superstring theories, the E 6 GUT gauge group occurs in four dimensions when we 
start with E 8 x E 8 in ten dimensions. The fermions are placed in a 2 7- dimensional represen- 
tation of E 6 . In such model, for each generation one would have new fermions including an 
isosinglet charge — | vector quark. 

In this article, we discuss the B meson radiative decay in the context of a generic vector 
quark model and show that the experimental data can be used to constrain the mixing 
angles. In vector quark models, due to the mixing of vector quarks with ordinary quarks, 
the Kobayashi-Maskawa (KM) matrix of the charged current interaction is not unitary. The 
internal flavor independent contributions in the W exchange penguin diagrams no longer 
cancel among the various internal up-type quarks. In addition, the mixing also generates 
non-zero tree level FCNC in the currents of Z boson and that of Higgs boson, which in turn 
gives rise to new penguin diagrams due to neutral boson exchanges. All these contributions 
will be carefully analyzed in this paper. Leading logarithmic (LL) QCD corrections are also 
included by using the effective Hamiltonian formalism. The paper is organized as follows: In 
section 2, we review the charged current interaction and the FCNC interactions in a generic 
vector quark model. Through the diagonalization of mass matrix, the non-unitarity of KM 
matrix and the magnitude of the FCNC can both be related to the mixing angles between 
vector and ordinary quarks. In section 3, various contributions to B meson radiative decays 
are discussed in the vector quark model. In section 4, we discuss constraints on the mixing 
angles from the new data on B radiative decays and from other FCNC effects. There are 
many previous analyses on the same issue. We shall make detailed comparison at appropriate 
points (mostly in section 3.) of our discussion. Most vector quark models in the literature 
are more complicated than the one we considered here. 

2 Vector Quark Model 

We consider the model in which the gauge structure of SM remains while one charge — | and 
one charge | isosinglet vector quarks are introduced. Denote the charge — | vector quark as 
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D and the charge | vector quark as U . Large Dirac masses of vector quarks, invariant under 
SU(2) L , naturally arise: 



M^UlUr + U R U L ) + M D {D L D R + D R D L ) 



(1) 



All the other Dirac masses can only arise from SU(2)l symmetry breaking effects. Assume 
that the weak SU(2) gauge symmetry breaking sector is an isodoublet scalar Higgs field 0, 
denoted as 



\ 
I 



= I „ I = I 1 

We can express the neutral field h in terms of real components: 

h° = H + i x . 



(2) 



(3) 



The conjugate of <fi is defined as 

1)0* 

-<fr J \ -(/>- 

Masses for ordinary quarks arise from gauge invariant Yukawa couplings: 

3 W.» 



V2 



/ 



(4) 



(5) 



In addition, gauge invariant Yukawa couplings between vector quarks and ordinary quarks 
are possible, which give rise to mixing between quarks of the same charge. For the model 
we are considering, these are: 



ff ^ L D R <P - ft ft L Uj - ft* 4> ] D R ^ L - ft* ftU R ^ L 



(6) 



In general, U will mix with the up-type quarks and D with down-type quarks. It is thus 
convenient to put mixing quarks into a four component column matrix: 



U L ,R 

t>L,R 
U T .„ 



d L ,R 
S L,R 

b L , R 
D L , R 



(7) 



where a — 1, 2, 3, 4. All the Dirac mass terms can then be collected into a matrix form: 

d' L M d d' R + d' R M\d' L and u' L M u u' R + u' R M[u' L . (8) 
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In this article, we use fields with prime to denote the weak eigenstates and those without 
prime to denote mass eigenstates. M. u ,d are 4x4 mass matrices. Since all the right- 
handed quarks, including vector quark, are isosinglet, we can use the right-handed chiral 
transformation to choose the right handed quark basis so that Ul, D L do not have Yukawa 
coupling to the ordinary right-handed quarks. In this basis, Aid and Ai u can be written as 



M D ) \ M v 

with 

M u,d = ^/M) J l,d = ~^fu,d ( 10 ) 

M^u (with hats) are the standard 3x3 mass matrices for ordinary quarks. J^ u is the 
three component column matrix which determines the mixings between ordinary and vector 
quarks. We assume that the bare masses M v D are much larger My/. With M UD factored 
out, Aid,u can be expressed in terms of small dimensionless parameters a, b: 





M d = M D \ , Ai u = M v . (11) 



The mixing matrix U^' d of the left-handed quarks and the corresponding one U^ d for 
right-handed quarks, defined as, 

U L,R = ^L,R U L,R) d L R = U L R d L R: (12) 

are the matrices that diagonalize M. u ^Ai\^ an d -M. ud JA U! d respectively. Hence the mass 
matrices can be expressed as 

M u = U u L m u U^ M d = U d L m d U d J (13) 

with m U) d the diagonalized mass matrices. The diagonalization can be carried out order 
by order in perturbation expansion with respect to small numbers a and b. For isosinglet 
vector quark model, the right-handed quark mixings are significantly smaller. The reason is 
that M^M d is composed of elements suppressed by two powers of a or b except for the (4, 4) 
element. As a result, the mixings of D R with dR,SR,bu are also suppressed by two powers 
of a or b. On the other hand, it can be shown that the mixings between Dl and bL,SL,di 
are only of first order in a or b. To get leading order results in the perturbation, one can 
assume that Ur = I. For convenience, write Ul as 



U L=\ *r ~ • ( 14 ) 
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where K is a 3 x 3 matrix and R, S are three component column matrices. To leading order 
in a and b, T is equal to 1. if equals the unitary matrix that diagonalizes aa) . The columns 
R and S, characterizing the mixing, are given by 

R = b, S = -Kb. (15) 

Now we can write down the various electroweak interactions in terms of mass eigenstates. 
The Z coupling to the left-handed mass eigenstates are given by 



9 r, fT » 



C z = —Z— Z M ( 4 - sin 2 e w J^), (16) 

COS fyj/ 



J£ = uOTtX + = ^l(^)7^l - ^d L (z d )^d L (17) 

The 4x4 matrices z are related to the mixing matrices by 

(z«) = C/^ a ^« 

(z rf ) = C/?a z l# (18) 

with a z = Diag(l, 1, 1, 0). Note that the matrix z is not diagonal. Flavor Changing Neutral 
Current (FCNC) is generated by the mixings between ordinary and vector quarksp], |3|, 
The charged current interaction is given by 

C w = j=(W-r + + W^Jn, (19) 
r- = u' L a w ^d' L = u L V>fd L (20) 

where a w = Diag(l, 1, 1, a) is composed of the Clebsch-Gordon coefficients of the correspond- 
ing quarks. For an isosinglet vector quark, a = 0. The 4x4 generalized KM matrix V is 
given by: 

V = Ut ] a w U d L . (21) 

The standard 3x3 KM matrix V KM is the the upper-left submatrix of V. Neither V nor V KM 
is unitary. Note that the non-unitarity of V is captured by two matrices 

(VV) = TjQa&J* 

(VV^) = U^a 2 w U u L . (22) 

In the model we are considering, these two matrices are identical to z u,d of the FCNC effects 
in Eq. [18] since is equal to a z . Indeed 



V*V = (z d ), V\rf = (z u ) (23) 
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This intimate relation between the non-unitarity of W charge current and the FCNC of Z 
boson is important for maintaining the gauge invariance of their combined contributions to 
any physical process. 

The off-diagonal elements of these matrices, characterizing the non-unitarity, is closely 
related to the mixing of ordinary and vector quarks. The off-diagonal elements are of order 
a 2 or b 2 . To calculate it, in principle, the next-to- leading order expansion of K, denoted as 
K 2 , is needed. In fact 

(V^ = (K d + Kf)^ + a 2 {b d )i{b d )* (24) 

Fortunately, by the unitarity of the mixing matrix U d , the combination k d + k d is equal 
to -(6 d )(6 d )t. 

k d 2 + kf = -(b d )(b d y (25) 

Thus the off-diagonal elements can be simplified 



(m, = (-l+« 2 )fe(M- 



(26) 



For isosinglet vector quark, a = 0. 

The Yukawa couplings between Higgs fields and quarks in weak eigenstate can be written 
in a matrix form as 



- -j= (^' L a z M d d' R (f) + d' R M d a z i/j' L ^ + tl)' L a z M u u'J> + u R Mta z ip' L ft 



(27) 



Note that a z is added to ensure that the left handed isosinglet vector quarks do not par- 
ticipate in the Yukawa couplings. The Yukawa interactions of quark mass eigenstates and 
unphysical charged Higgs fields ± are given by 



!J <u(m u VL - Vm d R)d] 0+ + 9 \d(-m d V ] L + V ] m u R)u] ^ 



while those of Higgs boson H and unphysical neutral Higgs field x by 

d(m d z d L + z d m d R)d + u{m u z u L + z u m u R)u\ H 
d(-m d z d L + z d m d R)d + u{m u z u L - z u m u R)u x° 



2M W 
~2M, 



w 



(28) 

(29) 
(30) 



3 B Meson Radiative Decay 

The B — > X s 7 decay, which already exists via one-loop VY-exchange diagram in SM, is known 
to be extremely sensitive to the structure of fundamental interactions at the electroweak scale 
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and serve as a good probe of new physics beyond SM because new interaction generically 
can also give rise to significant contribution at the one-loop level. 

The inclusive B — > X s ^ decay is especially interesting. In contrast to exclusive decay 
modes, it is theoretically clean in the sense that no specific low energy hadronic model is 
needed to describe the decays. As a result of the Heavy Quark Effective Theory (HQET), the 
inclusive B meson decay width T(B — > X s 7) can be well approximated by the corresponding 
b quark decay width T(b — > 57). The corrections to this approximation are suppressed 
by 1/t7t| H and is estimated to contribute well below 10% || [/J. This numerical limit is 
supposed to hold even for the recently discovered non-perturbative contributions which are 
suppressed by l/m 2 c instead of I /ml |5|. In the following, we focus on the dominant quark 
decay b — > 57. 

In SM, b — *■ 57 arises at the one loop level from the various W mediated penguin diagrams. 
The number of diagrams needed to be considered can be reduced by choosing the non-linear 
R,£ gauge as in [10|. In this gauge, the tri-linear coupling involving photon, W boson and the 
unphysical Higgs field + vanishes. Therefore only four diagrams , as in Fig. 1. contribute. 




Figure 1: Charged boson mediated penguin. 
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The on-shell Feynman amplitude can be written as 

iM(b -> si) = y^L± £ V ib V: s F 2 {x l )q^ v 7sa^{m b R + m s L)b (31) 

% 

with Xj = mj/Myy The sum is over the quarks u, c and £. The contributions to F 2 from the 
four diagrams are denoted as fYi , with the subscript 1 used to denote the diagrams with 
photon emitted from internal quark and 2 those with photon emitted from W boson. The 
functions /'s are given by 

ff{x) = e t {Ux)-Ui(x) + lUx)} , (32) 

f 2 w (x) = Z_ 1 ( x )-fa( x ) + 2t 1 {x)-iUx) , (33) 

ft(x) = l ei x[Ux)+Ux)} , (34) 

ft{x) = JxKo(x)-6(x)] • (35) 

Here the functions £(x) are defined as 

„i z n+1 dz lnx + (1 — x) + ■ ■ • + 

UX) = 7o l + (x-l)z = (1 - x)™+ 2 ^ ' (36) 

and 

s /' _^L_ = -J21. (37) 
JO 1 + (x — l)z 1 — X 

F 2 (x) is the sum of these functions and is given by 

= /rw + /rw + ftw + rfw = - 5r^ lni + 1 (38) 

For light quarks such as u and c, with x, — > 0, the first two terms on the right hand side 
are small and can be ignored, F 2 (x UtC ) is dominated by the x independent term ||. However 
these mass-independent terms get canceled among the up-type quarks due to the unitarity 
of KM matrix in SM 

E v*K = (39) 

% 

After the cancelation, the remaining contributions are essentially from penguins with internal 
t quark. 



It is convenient to discuss weak decays using the effective Hamiltonian formalism I 1 



12], which is crucial for incorporating the QCD corrections to be discussed later. In this 



formalism, the W and Z bosons are integrated out at matching boundary Myy- Their effects 
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are compensated by non-renormalizable effective Hamiltonian operators. The important 



dim-6 operators relevant for b — > 57 include 12 operators 



H* = -^KV*EC i (jM)O i 
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(40) 



Qi 
Q2 

Qa 

Q 5 

g 6 

Qi 
Qs 
Q9 
Q10 

g 7 

Qg 



(c a bf3)v-A(S/3C a )v-A , 
(c a 6a)y- J 4(S/3C^)y-A 

1 

(s a b/3)v-A ^2(qpqa)v-A 

q 

{s a b a )v~A ^{qpqp)v+A 

q 

(s a b/3)v-A ^2(q/3qa)v+A 



-(s a b a ) V -Aj2 e q(Wp)v+A 
A q 
3 

2' 

3 

2' 



q 

s a b a ) v - A e q(qpqp)v-A 



-{ s J>p)v-a e q(qpqa)v-A 

1 q 

s a a^[m b (l + 75 ) + m s {\ - -y 5 )]b a F, 



8tt 2 
9s 

8tt 2 



s a a^[m b {\ + 75) + m s (l - l5 )](T^)b a G A 



\1V 



(41) 



In Standard Model, the electroweak penguin operators Qr, . . . , Q10 are not necessary for a 
leading order calculation in 0(a). However, we will show later that in the vector quark 
model, FCNC effects exist as a linear combination of Q7, . . . , Q w . This effect will mix with 
Q 1 through RG evolution. 

The Wilson coefficients Cj at /i — M w are determined by the matching conditions when 
W, Z bosons and t quark are integrated out. To the zeroth order of a s and a, the only 
non-vanishing Wilson coefficients at /j, — M w for the above set are C 2 , C 7 , C G . C 2 is given 
by 

C 2 (M W ) = - V ; s V cb /V*V tb . (42) 



It is equal to one if the KM matrix is unitary and V* s V ub is ignored. C 7 at the scale M w is 
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given by the earlier penguin calculations 

C* M (M W ) = £ V ib V*F 2 (xi) = -\mxt) ~ -0.193 . (43) 

t.s tb i 

The numerical value is given when m t = 170 GeV. Here the function D' Q is defined as (jnj [T3| 

, » "I - Toe x (2 3^-1 , / , ,\ 

D » {x) s - 12(1 + l^r ( 44 > 

C 7 retains only the mass-dependent contribution from the penguin diagram with internal t 
quark. The mass-independent terms get cancelled, due to unitarity, among the three internal 
up type quarks. The mass-dependent contributions from the penguin diagrams with internal 
u and c quarks (they are small anyway) appear both in the high energy and low energy 
theories and get canceled in the matching procedure. Similarly, in SM the b — > sg transition 
arises from W exchange penguin diagrams which induce Qq. Since the gluons do not couple 
to W bosons, the gluonic W boson penguin consists only of two diagrams, which are given by 
/i with Q replaced by one. With the mass-independent contribution canceled, the Wilson 
coefficient Co can be written as 

C G M (M W ) = ~E' (x t ) ~ 0.096 , (45) 



The function E' is defined as |L3 



rn// \ x(x 2 -5x-2) 3 x 2 .. . 

4(1 -x)» + 2(T^r ]nx - (46) 

It is well known that short distance QCD correction is important for b — > S7 decay and 
actually enhances the decay rate by more than a factor of two. These QCD corrections can 
be attributed to logarithms of the form a"(m(,) \og m (rrib/Mw). The Leading Logarithmic 
Approximation (LLA) resums the LL series (m = n). Working to next-to-leading-log (NLL) 
means that we also resum all the terms of the form a s (mft)a™(mb) log 71 (rrif, / Mw) ■ The QCD 
corrections can be incorporated simply by running the renormalization scale from the match- 
ing scale fi = Myy down to m& and then calculate the Feynman amplitude at the scale m&. 
The anomalous dimensions for the RG running have been found to be scheme dependent 
even to LL order, depending on how 75 is defined in the dimensional regularization scheme. 

It has also been noticed [15j that the one- loop matrix elements of Q5, Qq for b — > S7 are 
also regularization scheme dependent. The matrix elements of Q^^ vanish in any four di- 
mensional regularization scheme and in the 't Hooft-Veltman (HV) scheme but are non-zero 
in the Naive dimension regularization (NDR) scheme. This dependence will exactly cancel 



= Octl3 
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the scheme-dependence in the anomalous dimensions and render a scheme-independent pre- 
diction. We refer to ]TT, T5\ for a review and details. In following, we choose to use the HV 



scheme. [] 

In the HV scheme, only Q 1 has a non-vanishing matrix element between b and 37, to 
leading order in a s (m b ), . Thus we only need C 7 (m b ) to calculate the LLA of b — > 37 decay 
width. For m t = 170 GeV, m b = 5 GeV, af\M z ) = 0.117 and in the HV scheme, C 7 (m 6 ) 
is related to the non-zero Wilson coefficients at My/ by [JTTJ, [L4|, |15 



C 7 (m b ) = 0.698 C^Mw) + 0.086 C G (M W ) - 0.156 C 2 (M W ). 
The b — > s-y amplitude is given by 

M (b - 87) = -V tb V t :^C 7 (m b )(Q 7 ) tiee (47) 

To avoid the uncertainty in m b , it is customary to calculate the ratio R between the radiative 
decay and the dominant semileptonic decay. The ratio R is given, to LLA, by [[14 



R = iff - S1 \ = jrh-^h \Vt* s V tb C 7 (m b )\ 2 . (48) 
Y[b^cev e ) \V cb \ TTg(z) 

Here the function g(z) of z = m c /m b is defined as 

g(z) = 1 - 8z 2 + 8z 6 - z 8 - 24z 4 In z (49) 

In the vector quark model, deviations from SM result come from various sources: (1) 
charged current KM matrix non-unitarity, (2) Flavor Changing Neutral Current (FCNC) 
effects in neutral boson mediated penguin diagrams, and (3) the W penguin with internal 
heavy U vector quark. Since the last one can be incorporated quite straight-forwardly, we do 
not elaborate on this contribution which will not be relevant for models without the U quark. 
We concentrate on the first two contributions, which have been discussed in Refs. |T6|, |T7| , |T8"| . 
Here we make a more careful and complete analysis which supplements or corrects these 



earlier analyses. Refs. |L6[ have calculated effects due to non-unitarity of the KM matrix and 
effects due to the Z mediated penguin in the Feynman gauge, however, their analysis did not 
include the FCNC contribution from the unphysical neutral Higgs boson, which is necessary 
for gauge invariance. The Higgs boson mediated penguins were also ignored. On the other 



1 It is customary in the literature to introduce ]l4j the scheme independent "effective coefficients" for 
Q 7 , Qg- These coefficients are defined so that <5 7 , Qg are the only operators with non-zero one loop matrix 
element for the process b — ► sj(g). The "effective coeffients" are hence identical to the original Wilson 
coefficients in HV scheme. 
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hand, Ref . [fTTf , while taking the unphysical Higgs boson into account, did not consider effects 
due to non-unitarity of the KM matrix, which gives the most important contribution. None 
of the above treatments, except Ref. 1T3] , included QCD corrections. 

Our strategy is first to integrate out the vector quark together with W and Z bosons at 
the scale My/. As shown above, the KM matrix is not unitary in the presence of an isosinglet 
vector quark. Therefore the mass-independent contributions in Eq. (|38]) from the magnetic 
penguin diagrams with various up-type quarks no longer cancel. This contribution is related 
to the short distance part of loop integration, i.e. when the loop momenta are large so that 
the quark mass which appears in the propagator can be ignored. In the formalism of the 
low energy effective Hamiltonian, it can be shown that these mass independent contributions 
never arise if the theory is renormalized using DRMS or similar schemes. By dimensional 
analysis, it is clear that the corresponding diagrams calculated using effective Hamiltonians 
are always proportional to the square of the loop quark mass. When we match the two 
calculations at the Mw scale, the mass-independent contributions must be compensated by 
new terms in Wilson coefficients. This is consistent with the notion that the effective field 
theory formalism separates the short distance physics encoded in Wilson coefficients from 
the long distance physics parameterized by matrix elements of the effective Hamiltonian. 
Such separation enables us to calculate effects of new physics by simply calculating Wilson 
coefficients perturbatively at the matching boundary. The matching results serve as initial 
conditions when Wilson coefficients run to a lower scale by renormalization group. Since the 
vector quarks have been integrated out, the anomalous dimensions are not affected by the 
new physics. 

Following this procedure, we calculate the extra contributions to the Wilson coefficient 
C 1 from non-unitarity: 

gnOgg 23 5 23 
V t * s V tb 36 V t * s V tb 36 • 1 ' 

The parameter S, one of the off-diagonal elements of the matrix V'V, characterizes the 
non-unitarity: 

S = (VW) 23 = z sh (51) 

The b — > 57 transitions also arise from FCNC Z boson and Higgs boson mediated penguin 
diagrams as in Fig. 2. The FCNC contribution to C 1 (Mw) can be denoted as follows: 



M + fs, b + O + =r£{fE + fh + fE) (52) 



VtbV *w •",<>• — VtbV * s 

For the sake of gauge invariance, f z needs to be considered together with f x . The Z boson 
penguins consist of internal charge — | quarks. The contribution from internal i = b,s quark, 
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fs,b, is S iven h J (Vi = ™ 2 /M 2 z ): 



fb = -\e d {(-\-e d sm 2 9 w )[2Uy b )-3^(y b )+Uyb)] 
+ e d sin 2 9 W [ 4£ (y b ) - 4& (y b )]} 

fs = -\e d {(-\-e d sm 2 9 w )[2Uys)-3^(y s )+Uys)] 

777 1 

m b J 



(53) 



(54) 



The last term in has a mass insertion in the s quark line. It is suppressed by m s /m b and 



will be ignored. The calculation is similar to that of f w . For a consistent approximation, 




a 



(b) 



Figure 2: Neutral boson mediated penguin diagrams. 



the two variables y b and y s , which are the ratios of m 2 , m 2 s to M§, are also set to zero. Hence 
fb+fs ~ -|e d {(-i- e(i sin 2 M[4eo(0)-66(0) + 26(0)] 



+e,sin 2 ^[4^o(0)-4ei(0)]} 
= -I-l sin 2 ^^-0.12 
The Z-mediated penguin diagram with internal D quark can also be calculated. 



fE = \e A [ 2Uvd) - 36 (y D ) + Hvd) ; 



(55) 
(56) 

(57) 



It approaches zero when yo — > oo and thus /p is negligible in large limit. For a gauge 
invariant result, the unphysical neutral Higgs x mediated penguin needs to be considered 
together with the Z boson penguin. In the non-linear Feynman gauge we have chosen, the 
mass of x is equal to M z . The calculation is very similar to the ^ penguin. For internal 
s, b, D quarks, the contributions f* b D are given by 



fi = ^yAiiivi) + &{yi) 



e d 

—z-yi 



v 2 - y t ) lny t - -yf + 4y 2 - — Vi + — 



Vi) 1 - 



(58) 
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It is obvious that the light quark contributions are suppressed by the light quark masses and 
thus negligible. The situation is quite different for the heavy D quark. As an approximation, 
for y D — > oo, /q — - > —^44 ~ —0.035. This contribution, comparable to the Z mediated 
penguin f z from light quarks, has been overlooked in previous calculations |16|. Since the 
quark D may not be much heavier than Z boson, we expand in powers oil/yn and keep 
also the next leading term. 



144 36 y D y z D 



The Higgs boson H mediated penguin is similar to that of unphysical Higgs \: 



e d 

Wi 

8 



7 15 
[—2 + 3wi) Inwi + -wf - <owf + —w. 



(60) 



where w% = mf/M^. Similar to the x penguin, ff b can be ignored since m s , nib <C tuh- For 
f§, we again expand it in powers of 1/wd and keep up to the next leading term: 



The leading term is +0.048, again comparable to the Z penguin. 

Put together, the Wilson coefficient C 1 (Mw) in the vector quark model is given by 

C y (M w ) = C™(M W ) + + -£-<J? + /* + ff + f b z + /* + f b H ) 

VtbVts 36 VtbVts 



VtbK 



(fE + M + fS) 



(MW) + V^ S U-g"^ 8111 ^ + 7^ + T44 - 36^-144 + T8^J 
^-0.193 + ^^x0.506 . (62) 

VaVZ 

Here we have used the unitarity relations z^bz\ s = — | 1 2 ^ s & ~ —z s b to leading order in 
FCNC due to the unitarity of Iff and 5 = z s b from Eq. (51). In the above numerical 
estimate we took y D , wd to infinity. 

Similarly the Wilson coefficient of the gluonic magnetic-penguin operator Qq is modified 
by the vector quark. In the vector quark model, the mass-independent term will give an extra 
contribution |5 if the KM matrix is non-unitary [|i~0] . The FCNC neutral boson mediated 
gluonic magnetic penguin diagrams are identical to those of the photonic magnetic penguin, 
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except for a trivial replacement of Qd by color factors, since photon and gluons do not couple 
to neutral bosons. Cq(M w ) in the vector quark model is given by 

C G (M W ) = C™(M W ) + -i-i - 3-%(/f + ft + /f + /* + /* + /f ) 

-3^ff(/I + /£ + /£) 

„c M , , , s z s b ( 1 1 1 . 9 . 5 5 11 7 1 1 \ 

= C^ M (M W ) + — — - + - + - s n 2 9 W + + 

G V tb V t * s \3 3 9 W 24y D 48 12 y D 48 6w D j 

-> -0.096 + x 0.733 . (63) 

The above deviation from SM does not include QCD evolution. We can incorporate LL 
QCD corrections to these deviations in the framework of effective Hamiltonian. The key is 
that the deviation from Standard Model is a short distance effect at the scale of M w and Mq. 
It can be separated into the Wilson coefficients at the matching scale, as we just did. The 
evolution of Wilson coefficients, which incorporates the LL QCD corrections, is not affected 
by the short distance physics of vector quark model and all the anomalous dimensions used 
in SM calculation still valid here. One only needs to use the corrected Wilson coefficients at 
fi = M w and in so doing we resum all the terms of the form z sb ai™(m;,) log™ (rrib/M w ). 

One subtlety in the vector quark model is that the quark mixing will generate Flavor 
Changing Neutral Current that couple to Z boson, which in turn gives rise to Z boson 
exchange interaction. This interaction is represented by an Effective Hamiltonian which is 
a linear combination of strong penguin operator Q 3 and electroweak penguin operator Qt^: 

G / 1 \ 

H NC = 2 ^=Zsb (--) (XMv-aXX^ - e q sin 2 9 w )(qpqp)v± A 

(64) 



v/2 S °V 2 
G F 



V2 Zs 



12 2 

~ o sin2 e wQi + o( X ~ sin2 °w)Qg 
6 3 3 



which gives additional non-zero Wilson coefficients: 

Z s b 1 



C 3 (M W ) 



VtbV t * s 6 ' 



C 7 (M W ) = -^^sin 2 ^, 

C 9 (M W ) = -^|(l-sin 2 M- (65) 

VtbVts 6 

To LL, the strong penguin and electroweak penguin operators could mixing among them- 
selves and also with Q 7 and Qg- This will generate an additional LL QCD correction to 
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b — > s 7 decay in the vector quark model. The crucial Wilson coefficient C^(m b ) obtain 
additional contributions: 

C y (m b ) = 0.698 C^Mw) + 0.086 C G {M W ) - 0.156 C 2 {M W ) 

+0.143 C 3 (M W ) + 0.101 C 7 (M W ) - 0.036 C 9 (M W ). (66) 

This FCNC LL QCD corrections is about one fifth the FCNC contribution of Z boson 
mediated penguin. The detail of the RG running calculation is given in the appendix. 

The correction to ratio R in the vector quark model, including its LL QCD corrections, 
is given by 

6ct 

AR = 1 \ it/ 12 x °- 307 x Re [KVtbZsb] = 0.23Rez s6 
to leading order in 8. In this result, the difference between V t * s V t b and —V* s V cbl i.e. 

V: s V cb = z sb - V*V tb7 (67) 

has been taken into account. We use the value | V^V^| 2 /| V cb \ 2 = 0.95. 



4 Constraints 

The inclusive B — > X s ^ branching ratio has been measured by CLEO with the branching 



ratio 19, 20 



B(B -> X s7 ) EXP = (3.15 ± 0.54) x 10" 4 (68) 

This branching ratio could be used to constrain the mixing in the vector quark model. We 
calculate the vector quark model deviation to leading logarithmic order, ie. all the terms 
of the form z s i ) a^{m b ) log" (m b /Mw). The Standard Model prediction to leading logarithmic 
order is |PT| 

B[B -»■ X s7 ) LO = (2.93 ± 0.67) x 10" 4 (69) 
The difference between the experimental data and the Standard Model LO prediction, 

B(B -> X s7 ) EX p - B{B -> X s7 ) NLO = (0.22 ± 0.86) x 10" 4 (70) 

It gives a range of possible vector quark model deviation and hence on z sb (with the input 
B(B -> X c eu) = 0.105): 

- 0.0027 < z sb < 0.0045 (71) 
^> The SM prediction up to next-to- leading order has been calculated in Ref.0, with 



=Octl3 
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the result 

B(B -> X s7 ) NLO = (3.28 ± 0.33) x (72) 

Ref. || later did a new analysis, which discards all corrections beyond NLO by expanding 
formulas like Eq.(f4"8"|) in powers of a s , and reported a slightly higher result: 

B(B -> X,7) NLO = (3.60 ± 0.33) x 10" 4 (73) 

Here we also use these new next-to-leading order SM calculations and the leading order vector 
quark model correction to constraint Z sb . To be consistent in the estimate of the theoretical 
errors, however, a full next-to-leading order calculation of the vector quark model matching 
correction is required. The difference between the experimental data and the Standard Model 
NLO prediction, with the errors added up directly, is 

B{B -> X s7 ) EX p - B(B -> X s7 ) NLO = (-0.13 ± 0.63) x 10- 4 [7] 

(-0.45 ± 0.63) x 10^ 4 @ (74) 

It gives a constraint on z sb . 

- 0.0032 < z sb < 0.0021 [7] 

-0.0045 < z sb < 0.0007 § (75) 

The previously strongest bound on z sb is from Z-mediated FCNC effect in the mode 
B^Xfi+fi- §: 

- 0.0012 < z sb < 0.0012 (76) 

Our new bound is as strong as that from FCNC. It shows that even though the vector quarks 
contribute to the radiative decay rate through one loop, as in SM, the data could still put 
strong bound. 

On the other hand, in models like Ref. [|]], operators of different chiralities such as 

°'i = -^s a a^[m b {l - 75 ) + m B (l + j 5 ))b a F^ , 
O'g = ^s a a^[m b (l - 75 ) + m.(l + l,)](T*)b a G% , (77) 

occurs via the new interaction. Our study can be extended to these models too. However, 
the new amplitude for b — > S7 belongs to a different helicity configuration in the final state 
and it will not interfere with the SM contribution. Consequently, the constraint obtained 
from b — > s 7 in these models is less stringent than that from B — > Xfi + {i~. 
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In the upcoming years, much more precise measurements are expected from the upgraded 
CLEO detector, as well as from the 7?-factories presently under construction at SLAC and 
KEK. The new experimental result will certainly give us clearer evidence whether the vector 
quark model is viable. 

DC wishes to thank T. Morozumi and E. Ma for discussions. WYK is supported by a 
grant from DOE of USA. DC and CHC are supported by grants from NSC of ROC. CHC's 
support by NSC is under contract No.NSC 88-2112-M-003-017. 



Appendix 



The RG equations for the Wilson coefficients C r = (Ci(//), 
written as 



Cio(aO) an d C7; Cg can be 



d 



din fx 

d 



din// 

d 

din // 



cm 

cm 



a 



^iJcm 



Ots_ 

Air 



(Ay • cm + i^cm + tg 7 Cg(/^)) 

((3 g -CM+1ggC g (h)) 



(7* 



The 10 by 10 submatrix 7 r can be found in ||15|| . The anomalous dimension matrix entries 
0L~o° are extracted from the results of (3 z ~q in Ref. In the HV scheme, (3 lG are 

given by: 





4 N 2 -l 2e u N 2 -l 



\ 



27 2N e d 2N 

116 N 2 -l 

27 2N 
4/ N 2 -l n JN 2 -l 



2fi 



2N 



27 2N 

8 iV 2 -l 

3 2N 
4/ jV 2 -l 1 o f N 2 -l 
"27 2N J 2N 

27V _ 



2JV 

58 Af 2 



-1 



9 

\ 9 Crf 27V 



2/V 



3/ 



2iV 



iV 2 -l 
2JV 



(3 



G 



\ 



UN 
9 



29 
9JV 



6 + 



29/ 
9JV 

-3/ 

I 25/ 

9JV 



9 9iV 

lijyi 

9 



-3/Ve d + £e d - fe d 



N 

-6e d - ^le d + fie, 

o 6, n«6{i "I - ~TT&d 



3 " a 3/V 
11 AT/ 



V 9e d + 



2_ 



(79) 



Here u and d are the numbers of active up-type quarks and down type quarks respectively, 
/ is the total number of active quark flavor. Between the scales m& and Mw, u = 2, d = 3, 



/ 



(e d d + e u u)/e d 



-1, / 



[e\d + e2 M )/e d = -11/3. For 377(3) color, TV = 3 . 
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